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1 Introduction 



Eugene Paul Wigner received the 1963 Noble prize in physics for his contri- 
butions to symmetry problems in physics. Even before the formulation of the 
present form of quantum mechanics, Wigner perceived the importance of sym- 
metry problems in quantum mechanics. In his book which was published in 
1931 PP, Wigner formulated the quantum theory of angular momentum in terms 
of the three-dimensional rotational symmetries. In so doing, he introduced group 
theoretical methods to physics. 

In 1932 |2], he published a paper concerning thermodynamic corrections to 
equilibrium systems, and introduced a phase-space picture of quantum mechan- 
ics. The phase-space distribution function he introduced in this paper is called 
the Wigner function. The Wigner function is an indispensable theoretical tool 
in quantum optics and in the foundations of quantum mechanics 

Frederick Seitz was Wigner's first student at Princeton. In their papers pub- 
lished in 1933 and 1934 4 , they worked out the symmetry of sodium crystals, 
in terms of group theory. They initiated application of quantum mechanics to 
matter. This field is now known as condensed matter physics. 

By 1936, the spin and isospin symmetries for nucleons were well established 
in physics. Wigner was interested whether those two separate symmetries come 
from one big symmetry. In so doing, he formulated the concept of supermul- 
tiplets Wigner's supcrmultiplet theory was the grandfather of the quark 
model which Gell-Mann formulated in 1963 [2]. 

In 1939, Wigner published his most fundamental paper dealing with space- 
time symmetries of relativistic particles [7]. In this paper, Wigner introduced 
the Lorentz group to physics. Furthermore, by introducing his "little groups," 
Wigner provided the framework for studying the internal space-time symme- 
tries of relativistic particles. Since this paper was so ahead of his time, it was 
rejected by three different journals before John von Neumann, then the editor 
of the Annals of Melodramatics, decided to publish it in his journal. The sci- 
entific contents of this paper have not yet been fully recognized by the physics 
community. We are writing this report as a continuation of the work Wigner 
initiated in this history-making paper. 

While particle physicists are still struggling to understand internal space- 
time symmetries of elementary particles, Wigner's Lorentz group is becoming 
useful to many other branches of physics. Among them is optical sciences, both 
quantum and classical. In quantum optics, the coherent and squeezed states 
are representations of the Lorentz group [3]. Recently, the Lorentz group is 
becoming the fundamental language for classical ray optics. It is gratifying to 
note that optical components, such as lenses, polarizers, interferometers, lasers, 
and multi-layers can all be formulated in terms of the Lorentz group which 
Wigner formulated in his 1939 paper. Classical ray optics is of course a very 
old subject, but we cannot do new physics without measurements using optical 
instruments. Indeed, classical ray optics constitutes Wigner's frontier in physics. 

The word "group theory" sounds like an abstract mathematics, but it is 
gratifying that Wigner's little groups can be formulated in terms of two-by-two 
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Table 1: Further contents of Einstein's E = mc 2 . Massive and massless par- 
ticles have different energy-momentum relations. Einstein's special relativity 
gives one relation for both. Wigner's little group unifies the internal space-time 
symmetries for massive and massless particles. The quark model and the parton 
model can also be combined into one covariant model. 
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matrices, while classical ray optics is largely a physics of two-by-two matrices. 
The mathematical correspondence is straight-forward. 

In order to see this point clearly, let us start with the following classic ex- 
ample. The second-order differential equation 

A^ + B d A +Cq{ t)=Fco s{ .t). (1) 

is applicable to a driven harmonic oscillator with dissipation. This can also 
be used for studying an electronic circuits consisting of inductance, resistance, 
capacitance, and an alternator. Thus, it is possible to study the oscillator 
system using the electronic circuit. Likewise, an algebra of two-by-two matrices 
can serve as the scientific language for two more different branches of physics. 

There are many physical systems which can be formulated in terms of two- 
by-two matrices. If we restrict that their determinant be one, there is a well 
established mathematical descipline called the group theory of SL(2,C). This 
aspect was noted in the study of Lorcntz transformations. In group theoretical 
terminology, the group SL(2, C) is the universal covering group of the group 
of the Lorentz group. In practical terms, to each two-by-two matrix, there 
corresponds one four-by-four matrix which performs a Lorentz transformation 
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in the four-dimensional Minkowskian space. Thus, if a physical system can be 
explained in terms of two- by- two matrices, it can be explained with the language 
of Lorentz transformations. Conversely, the system can serve as an analogue 
computer for Lorentz transformations. 

Optical filters, polarizers, and interferometers deal with two independent 
optical rays. They superpose the beams, change the relative phase shift, and 
change relative amplitudes. The basic language here is called the Jones matrix 
formalism, consisting of the two-by-two matrix representation of the SL(2,C) 
group IH] . The four- by- four Mueller matrices are derivable from the two-by- 
two matrices of SL(2,C). 

Para-axial lens optics can also be formulated in terms of two-by-two matrices, 
applicable to the two-component vector space consisting of the distance from 
the optical axis and the slope with respect to the axis. The lens and translation 
matrices are triangular, but they are basically representations of the Sp(2) group 
which is the real subgroup of the group SL(2, C) [TnUTT] . 

Laser optics is basically multi-lens lens optics. However, the problem here is 
how to get simple mathematical expression for the system of a large number of 
the same lens separately by equal distance. Here again, group theory simplifies 
calculations [T2). 

In multi-layer optics, we deal with two optical rays moving in opposite di- 
rections. The standard language in this case is the S-matrix formalism 
This is also a two-by-two matrix formalism. As in the case of laser cavities, the 
problem is the multiplication of a large number of matrix chains [141 115j . 

It is shown in this report that the two-by-two representation of the six- 
parameter Lorentz group is the underlying common scientific language for all of 
the instruments mentioned above. While the abstract group theoretical ideas 
make two-by-two matrix calculations more systematic and transparent in optics, 
optical instruments can act as analogue computers for Lorentz transformations 
in special relativity. It is gratifying to note that special relativity and ray optics 
can be formulated as the physics of two-by-two matrices. 

In Sec. |21 we discuss the historical significance of Winger's 1939 paper [7] 
on the Lorentz group and its application to the internal space-time symmetries 
of relativistic particles. In Sec. |31 we present the basic building blocks for 
the two-by-two representation of the Lorentz group in terms of the matrices 
commonly seen in ray optics. In Sees. 01 El 13 and [HI we discuss polarization 
optics, one-lens system, multi-lens system, laser cavities, and multi-layer optics, 
respectively. 

Since the Lorentz group is relatively new to many who study optics, we 
explain how it is possible to represent the group using two-by-two matrices in 
Appendix ^ we explain how the Lorentz group can be formulated in terms 
of four- by- four matrices. It is shown that the group can have six independent 
parameters. In Appendix [Bj we explain how it is possible to formulate the 
Lorentz group in terms of two- by- two matrices. It is shown that the four- by- four 
transformation matrices can be constructed from those two-by-two matrices. In 
Appendix^! it is noted that the four-by-four matrices are real, their two-by-two 
counterparts are complex. However, there is a three-parameter subgroup called 
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Sp(2). It is shown that the complex subgroup SU(1, 1) is equivalent to Sp(2) 
through conjugate transformation. 

The purpose of these appendixes is to give an introduction to group the- 
oretical methods used in this report and in the recent optics literature. In 
Appendix IdI it is shown that the Lie group method, in terms of the generators, 
is not the only method in constructing group representations. For the rotation 
group and the three-parameter subgroups of the Lorentz group, it is simpler to 
start with the minimum number of starter matrices. For instance, while there 
are three generators for the rotation group in the Lie approach, we can con- 
struct the most general form of the rotation matrix from rotations around two 
directions, as Goldstein constructed the Euler angles 

As for the references, we have not made attempts to list all the papers 
relevant to the present report. However the references are given in the papers 
in the refereed journals. 

2 Wigner's Little Groups 

If the momentum of a particle is much smaller than its rest-mass energy, the 
energy-momentum relation is E = p 2 /2mo + moc 2 . If the momentum is much 
larger than its rest-mass energy, the relation is E = cp. These two different 
relations can be combined into one covariant formula E 2 = m^c 4 +p 2 c 2 . This 
aspect of Einstein's E — mc 2 is also well known. 

In addition, particles have internal space-time variables. Massive particles 
have spins while massless particles have their helicities and gauge variables. 
Our first question is whether this aspect of space-time variables can be unified 
into one covariant concept. The answer to this question is Yes. Wigner's little 
group does the job. In addition, particles can have space-time extensions. For 
instance, in the quark model, hadrons are bound states of quarks. However, 
the hadron appears as a collection of partons when they it moves with speed 
close to the velocity of light. Quarks and partons seem to have quite distinct 
properties. Are they different manifestation of a single covariant entity? This 
is one of the main issues in high-enrgy particles physics. 

By "further contents" of Einstein's E = mc 2 , we mean that the internal 
space-time structures of massive and massless particles can be unified into one 
covariant package, as E 2 — TOqC 4 + p 2 c 2 does for the energy-momentum rela- 
tion. The mathematical framework of this program was developed by Eugene 
Wigner in 1939 7 . He constructed maximal subgroups of the Lorentz group 
whose transformations will leave the four-momentum of a given particle invari- 
ant. These groups are known as Wigner's little groups. Thus, the transforma- 
tions of the little groups change the internal space-time variables the particle. 
The little group is a covariant entity and takes different forms for the particles 
moving with different speeds. 

The space-time symmetry of relativistic particles is dictated by the Poincare 
group |7|. The Poincare group is the group of inhomogeneous Lorentz trans- 
formations, namely Lorentz transformations preceded or followed by space-time 
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translations. Thus, the Poincare group is a semi-direct product of the Lorentz 
and translation groups. The two Casimir operators of this group correspond to 
the (mass) 2 and (spin) 2 of a given particle. Indeed, the particle mass and its 
spin magnitude are Lorentz- invariant quantities. 

The question then is how to construct the representations of the Lorentz 
group which are relevant to physics. For this purpose, Wigner in 1939 stud- 
ied the maximal subgroups of the Lorentz group whose transformations leave 
the four- momentum of a given free particle [7]. These subgroups are called 
the little groups. Since the little group leaves the four-momentum invariant, 
it governs the internal space-time symmetries of relativistic particles. Wigner 
shows in his paper that the internal space-time symmetries of massive and mass- 
less particles are dictated by the little groups which are locally isomorphic to 
the three-dimensional rotation group and the two-dimensional Euclidean groups 
respectively. 

The group of Lorentz transformations consists of three boosts and three 
rotations. The rotations therefore constitute a subgroup of the Lorentz group. 
If a massive particle is at rest, its four-momentum is invariant under rotations. 
Thus the little group for a massive particle at rest is the three-dimensional 
rotation group. Then what is affected by the rotation? The answer to this 
question is very simple. The particle in general has its spin. The spin orientation 
is going to be affected by the rotation! If we use the four-vector coordinate 
(t, z, x, y), the four-momentum vector for the particle at rest is (tooc 2 , 0, 0, 0), 
and the three-dimensional rotation group leaves this four-momentum invariant. 
This little group is generated by 

/0 0\ /0 0\ /0 

[0 Oil -i 1 [000 

1 ~ 0' 2 ~ i oh 3 _ -i 
\0 -i 0/ \0 0/ \0 i 

(2) 

These are essentially the generators of the three-dimensional rotation group. 
They satisfy the commutation relations: 

[Ji, Jj] = ieijkJk- (3) 

If the rest-particle is boosted along the z direction, it will pick up a non-zero 
momentum component along the same direction. The above generators will also 
be boosted. The boost will take the form of conjugation by the boost matrix 

(cosh 77 sinh 77 \ 

sinhry cosh r] 1 

10' [ ' 

1/ 

This boost will not change the commutation relations of Eq.© for 0(3), and 
the boosted little group will still leave the boosted four- momentum invariant. 
Thus, the little group of a moving massive particle is still 0(3)-like. 
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It is not possible to bring a massless particle to its rest frame. In his 1939 
paper |7], Wigner observed that the little group for a massless particle moving 
along the z axis is generated by the rotation generator around the z axis, namely 
J3 of Eq.@, and two other generators which take the form 








i 











i 





i 


—i 








Vo 












No 



/0 < 

1 



\i -i 0, 



(5) 



If we use Ki for the boost generator along the i-th axis, these matrices can be 
written as 



with 



Ni = K\ - J 2 

/0 i 0^ 



1 
\0 0, 



N 2 = K 2 + Ji, 



Ki 



/0 V 





\i 0, 



(6) 



(7) 



The generators J3, N\ and N 2 satisfy the following set of commutation relations. 

[7Vi,iV 2 ] = 0, [J 3 ,Nt]=iN 2 , [J 3 ,N 2 ] = -iN 1 . (8) 

In order to understand the mathematical basis of the above commutation 
relations, let us consider transformations on a two-dimensional plane with the 
xy coordinate system. We can then make rotations around the origin and trans- 
lations along the x and y directions. If we write these generators as L, P x and 
P y respectively, they satisfy the commutation relations 16 



[P x ,Py] - 0, [L, P X ] = tPy, [L, Py] = -iP 3 



(9) 



This is a closed set of commutation relations for the generators of the E(2) 
group. If we replace N± and N 2 of Eq.JHJ) by P x and P y , and J 3 by L, the 
commutations relations for the generators of the £ , (2)-like little group becomes 
those for the 2£(2)-like little group. This is precisely why we say that the little 
group for massless particles are like E(2). 

It is not difficult to associate the rotation generator J3 with the helicity 
degree of freedom of the massless particle. Then what physical variable is as- 
sociated with the N\ and N 2 generators? Indeed, Wigner was the one who 
discovered the existence of these generators, but did not give any physical inter- 
pretation to these translation- like generators in his original paper 0. For this 
reason, for many years, only those representations with the zero-eigenvalues of 
the N operators were thought to be physically meaningful representations |17| . 
It was not until 1971 when Janner and Janssen reported that the transforma- 
tions generated by these operators are gauge transformations |18II19| . The role 
of this translation-like transformation has also been studied for spin- 1/2 par- 
ticles, and it was concluded that the polarization of neutrinos is due to gauge 
invariance |20l 121 j . 
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The 0(3)-like little group remains 0(3)-like when the particle is Lorentz- 
boosted. Then, what happens when the particle speed becomes the speed of 
light? The energy-momentum relation E 2 = m^c 4 + p 2 c 2 become E = pc. Is 
there then a limiting case of the 0(3)-like little group? Since those little groups 
are like the three-dimensional rotation group and the two-dimensional Euclidean 
group respectively, we are first interested in whether E(2) can be obtained from 
0(3). This will then give a clue to obtain the _E(2)-like little group as a limiting 
case of 0(3)-like little group. With this point in mind, let us look into this 
geometrical problem. 

In 1953, Inonu and Wigner formulated this problem as the contraction of 
0(3) to E{2) [22]. Let us see what they did. We always associate the three- 
dimensional rotation group with a spherical surface. Let us consider a circular 
area of radius 1 kilometer centered on the north pole of the earth. Since the 
radius of the earth is more than 6,450 times longer, the circular region appears 
flat. Thus, within this region, we use the E(2) symmetry group. The validity 
of this approximation depends on the ratio of the two radii. 

How about then the little groups which are isomorphic to 0(3) and E(2)l 
It is reasonable to expect that the i?(2)-like little group can be obtained as a 
limiting case for of the 0(3)-like little group for massless particles. In 1981, it 
was observed by Bacry and Chang [23] and by Ferrara and Savoy [23 that this 
limiting process is the Lorentz boost to infinite-momentum frame. 

In 1983, it was noted by Han et al that the large-radius limit in the the 
contraction of 0(3) to E(2) corresponds to the infinite-momentum limit for 
the case of 0(3)-like little group to i?(2)-like little group. They showed that 
transverse rotation generators become the generators of gauge transformations 
in the limit of infinite momentum |25| . 

Let us see how this happens. The J3 operator of Eq.J2|, which generates 
rotations around the z axis, is not affected by the boost conjugation of the B 
matrix of Eq.Q. On the other hand, the J\ and J2 matrices become 

N\ = lim e- 7, B- 1 J 2 B, N 2 = lim -e^B' 1 J X B, (10) 

77 — >oo 77 — »oo 

and they become Ni and N2 given in Eq. JSJ) . The generators Ni and N% are the 
contracted J2 and J\ respectively in the infinite-momentum. In 1987, Kim and 
Wigner studied this problem in more detail and showed that the little group 
for massless particles is the cylindrical group which is isomorphic to the E(2) 
group [2Ej- 

This completes the second row in Tabled where Wigner 's little group uni- 
fies the internal space-time symmetries of massive and massless particles. The 
transverse components of the rotation generators become generators of gauge 
transformations in the infinite- momentum limit. 

Let us go back to Table I given in Sec.^ As for the third row for relativistic 
extended particles, the most efficient approach is to construct representations 
of the little groups using the wave functions which can be Lorentz-boosted. 
This means that we have to construct wave functions which are consistent with 
all known rules of quantum mechanics and special relativity. It is possible 
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to construct harmonic oscillator wave functions which satisfy these conditions. 
We can then take the low-speed and high-speed limits of the covariant harmonic 
oscillator wave functions for the quark model and the parton model respectively. 
This aspect was extensively discussed in the literature |16|. and is beyond the 
scope of the present report. 



3 Formulation of the Problem 

Let us consider two optical beams propagating along the z axis. We are then 
led to the column vector: 



A\ exp (— i (kz — cut + 0i)) 
A2 exp (— i (kz — Lot + ^2)) 



(11) 



We can then achieve a phase shift between the beams by applying the two-by- 
two matrix: 

V * /2 ° ^ (12) 
If we are interested in mixing up the two beams, we can apply 



cos(0/2) -sin(0/2)\ 
sin(0/2) cos(0/2) J 



(13) 



to the column vector. 

If the amplitudes become changed by either by attenuation or reflection, we 
can use the matrix 

'e"/ 2 
e""/ 2 



(14) 



for the change. In this paper, we are dealing only with the relative amplitudes, 
or the ratio of the amplitudes. 

Repeated applications of these matrices lead to the form 



a (3 
7 8 



(15) 



where the elements are in general complex numbers. The determinant of this 
matrix is one. Thus, the matrix can have six independent parameters. 

Indeed, this matrix is the most general form of the matrices in the SL(2, c) 
group, which is known to be the universal covering group for the six-parameter 
Lorentz group. This means that, to each two-by-two matrix of SL{2, c), there 
corresponds one four-by-four matrix of the group of Lorentz transformations ap- 
plicable to the four-dimensional Minkowski space • It is possible to construct 
explicitly the four-by-four Lorentz transformation matrix from the parameters 
a, (3, 7, and 8. This expression is available in the literature ^Hjj an d we consider 
here only special cases. 
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As is shown in Appendix^] the Lorentz group includes Lorentz boosts along 
three different directions and rotations around those three different directions. 
A rotation around the z axis in the Minowskian space can be written as 



/l 

10 

cos (f) — sin ( 

\0 s\n(h cost 



(16) 



This four- by- four matrix corresponds to the two- by- two matrix of Ea. (|13[) . We 
use in this paper the four-vector convention (t, z, x, y). 
The rotation around the y axis can be written as 

/l 0\ 

cos 8 —sin 8 1 , 

sinfl cos(9 ' ^ ' 

Vo 1/ 

The two-by-two matrix of Eq. i|14|) corresponds to the Lorentz boost matrix 

(18) 



The four-by-four representation of the Lorentz group can be constructed 
from the two-by-two representation of the SL(2, c) group, which is known as 
the universal covering group of the six-parameter Lorentz group. This aspect is 
discussed in Appendix FBI 

Let us go back to Ea. l|15|l . the SL(2,c) group represented by this matrix 
has many interesting subgroups. If the matrices are to be Hermitian, then the 
subgroup is SU(2) corresponding three-dimensional rotation group. If all the 
elements are real numbers, the group becomes the three-parameter Sp(2) group. 
This subgroup is equivalent to SU(1, 1) which is the primary scientific language 
for squeezed states of light [23 Ej • 

We can also consider the matrix of Eq.JTSJ when one of its off-diagonal 
elements vanishes. Then, it takes the form 




exp (i(f>/2) 

7 exp (—i0/2) 



(19) 



where 7 is a complex number with two real parameters. In 1939 [7j, Wigner 
observed this form as one of the subgroups of the Lorentz group. He observed 
further that this group is isomorphic to the two-dimensional Euclidean group, 
and that its four-by-four equivalent can explain the internal space-time symme- 
tries of massless particles including the photons. 

In ray optics, we often have to deal with this type of triangluar matrices, 
particularly in lens optics and stability problems in laser and multilayer optics. 
In the language of mathemtics, dealing with this form is called the Iwasawa 
decomposition |28j . This aspect of the Lorentz group is dicussed in Appendix ICl 
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4 Polarization Optics and Interferometers 



In polarization physics, we are studying optical rays with two independent com- 
ponents. Thus the two-by-two matrix given in Sec ^ is directly applicable. We 
are quite familiar with Pauli spinors and Pauli matrices with three independent 
parameters. They deal with rotations. If we add Lorentz boosts, the group be- 
comes the six-parameter SL(2, c) group. The representation is still two-by-two. 

From the group theoretical point of view, the Jones matrices constitute the 
SL(2, c) group or the universal covering group for the Lorentz group. The 
Jones vectors and Jones matrices are nothing more and less than the S'L(2,C) 
representation and the SL{2) spinors. 

There is a specific procedure to construct the Minkowskian four-vectors and 
the four-by-four Lorentz transformation matrices. Then they become Stokes 
parameters and the Mueller matrices. 

If the Jones matrix contains all the parameters for the polarized light beam, 
why do we need the mathematics in the four-dimensional space? The answer 
to this question is well known. In addition to the basic parameter given by the 
Jones vector, the Stokes parameters give the degree of coherence between the 
two rays. 

Let us write Eq. as a Jones spinor of the form 

($::!)■ 

and construct the quantities: 

Sii =< V'lV'i. >j Su =< >, 

S 21 = < >> S 22 =< Tp* 2 4>2 > • (21) 

Then the Stokes vector consists of 

So = Sn + S22, Si = Sn — S22, 



S 2 = S12 + S21, S 3 = -t(Sn + S22). (22) 

The four-component vector (So, Si, S2, S3) transforms like the space-time four- 
vector (t, z, x, y) under Lorentz transformations. The Mueller matrix is therefore 
like the Lorentz-transformation matrix. 

As in the case of special relativity, let us consider the quantity 

M 2 = Si - Si - Si - Si. (23) 

Then M is like the mass of the particle while the Stokes four-vector is like the 
four- momentum. 

If M = 0, the two-beams are in a purely state. As M increases, the system 
becomes mixed, and the entropy increases. If it reaches the value of Sq, the 
system becomes completely random. It is gratifying to note that this mechanism 
can be formulated in terms of the four-momentum in particle physics [3]. 
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5 One-lens System 



In analyzing optical rays in para-axial lens optics, we start with the lens matrix: 

-i 1 // !)• (24) 

and the translation matrix 

T=(j *), (25) 

assuming that the beam is propagating along the Z direction. 
Then the one-lens system consists of 



1 2 2 \f 1 1 Zi 

1 J I -1// 10 1 



(26) 



(27) 



If we perform the matrix multiplication, 

I-Z2// z i + z 2 - ZlZ 2 /f 
-l/f l-zi/f 

If we assert that the upper-right element be zero, then 

zi z 2 f 

and the image is focussed, where z\ and z-i are the distance between the lens 
and object and between the lens and image respectively. They are in general 
different, but we shall assume for simplicity that they are the same: z\ = z 2 = z. 
We are doing this because this simplicity does not destroy the main point of our 
discussion, and because the case with two different values has been dealt with 
in the literature |15|. Under this assumption, we are left with 

i-z/f 2z- z y f \ 

-l/f 1-z/f J' [M) 

The diagonal elements of this matrix are dimensionless. In order to make the 
off-diagonal elements dimensionless, we write this matrix as 



sTz \(l-z/f zlf-2\(yfi 

l/y/z) I Zff l-z/f) \ \/y/z 



(30) 



Indeed, the matrix in the middle contains dimensionless elements. The negative 
sign in front is purely for convenience. We are then led to study the core matrix 

c=( x ~ i x ~?y (31) 

V X x — 1 I 
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Here, the important point is that the above matrices can be written in terms 
of transformations in the Lorentz group. In the two-by-two matrix representa- 
tion, the Lorentz boost along the z direction takes the form 

and the rotation around the y axis can be written as 

R(e) = ( co fJ / l\ 03) 

v V sm(0/2) cos(</>/2) / x ' 

and the boost along the x axis takes the form 

Y(v\ - ( cosh (x/ 2 ) sinh(x/2) \ , - 

[X) - \ sinh( X /2) cosh( X /2) J ' 1 ' 

Then the core matrix of Eq. I|31[) can be written as 

B{r,)R{cj>)B{-r,), (35) 



(36) 



or 

cos(0/2) -e-''sin(0/2)\ 
;"sin(0/2) cos(</>/2) )' 

if 1 < x < 2. If x is greater than 2, the upper-right element of the core is 
positive and it can take the form 

B{n)X{ X )B{-n), (37) 



or 

/ cosh( X /2) e-"sinh(x/2)\ 

l v e+"sinh(x/2) cosh(x/2) )' y ' 

The expressions of Eq. I|35|l and Eq. 137|l are a Lorentz boosted rotation and a 
Lorentz-boosted boost matrix along the x direction respectively. These expres- 
sions play the key role in understanding Wigner's little groups for relativistic 
particles. 

Let us look at their explicit matrix representations given in Ea. (|36|l and 
Eci. l38|l . The transition from Ea. (|36|l to Ea. l|38|l requires the upper right ele- 
ment going through zero. This can only be achieved through r\ going to infinity. 
If we like to keep the lower-left element finite during this process, the angle <f> 
and the boost parameter x have to approach zero. The process of approach- 
ing the vanishing upper-right element is necessarily a singular transformation. 
This aspect plays the key role in unifying the internal space-time symmetries 
of massive and massless particles. This is like Einstein's E = \J (pc) 2 + mf ) c A 
becoming E = pc in the limit of large momentum. 

On the other hand, the core matrix of Ea. H31(l is an analytic function of 
the variable x. Thus, the lens matrix allows a parametrization which allows 
the transition from massive particle to massless particle analytically. The lens 
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optics indeed serves as the analogue computer for this important transition in 
particle physics. 

^From the mathematical point of view, Eg. ((36(1 and Eg. 1(38(1 represent cir- 
cular and hyperbolic geometries, respectively. The transition from one to the 
other is not a trivial mathematical procedure. It requires a further investigation. 

Let us go back to the core matrix of Eo. ((31(l . The x parameter does not 
appear to be a parameter of Lorentz transformations. However, the matrix can 
be written in terms of another set of Lorentz transformations. This aspect has 
been discussed in the literature [TT] . 

6 Multi-lens Problem 

Let us consider a co-axial system of an arbitrary number of lens. Their focal 
lengths are not necessarily the same, nor are their separations. We are then 
led to consider an arbitrary number of the lens matrix given in Eq.{21|| and an 
arbitrary number of translation matrix of Ea. (|25|) . They are multiplied like 

Ti L\ T 2 L 2 T 3 L 3 Tjv Ljy, (39) 

where N is the number of lenses. 

The easiest way to tackle this problem in to use the Lie-algebra approach. 
Let us start with the generators of the Sp(2) group: 

Since the generators are pure imaginary, the transformation matrices are real. 

On the other hand, the L and T matrices of Eg . I(24|) and Eg. ((2511 are gener- 
ated by 

'.-(SS). *-(!!!)■ <«> 

If we introduce the third matrix 

Mo -")■ m 

all three matrices form a closed set of commutation relations: 

[X u X 2 ]=iX 3 , [X u X 3 ] = -iX u [X 2 ,X 3 ]=iX 2 . (43) 

Thus, these generators also form a closed set of Lie algebra generating real two- 
by-two matrices. What group would this generate? The answer has to be Sp(2). 
The truth is that the three generators given in Eq. 1(43(1 can be written as linear 
combinations of the generators of the Sp(2) group given in Eq.^UJ) JTOj. Thus, 
the Xi matrices given above can also act as the generators of the Sp(2) group, 
and the lens-system matrix given in Eq. I|39() is a three-parameter matrix of the 
form of Ea. HI 5(1 with real elements. 



14 



The resulting real matrix is written as 

c I) • <«> 

and is called the ABCD matrix. According to Bargamann [3U|, this three- 
parameter matrix can be decomposed into 



cos(a/2) -sin(a/2)\ ( e 1 ' 2 \ /cos(/3/2) -sin(/?/2)\ 
sin(a/2) cos(a/2) )\ e"^ 2 J ^ sin(/3/2) cos(/?/2) J' 



(45) 



which can be written as the product of one symmetric matrix resulting from 

(46) 

and one rotation matrix: 

(47) 



cos(a/2) -sin(a/2)\ ( e 1 ' 2 \ { cos(a/2) — sin(a/2) 
sin(a/2) cos(a/2) J \ e^/ 2 J I sin(a/2) cos(a/2) 



cos[(/3 - a) /2] - sin[(/3 - a) /2] 
sin[(/3- a)/2] cos[(/3 - a)/2] 

We can then decompose each of these two matrices into the lens and translation 
matrices. The net result is that we do not need more than three lenses to de- 
scribe the lens system consisting of an arbitrary number of lenses. The detailed 
calculations are given in Ref. |10j . 

7 Laser Cavities 

In a laser cavity, the optical ray makes round trips between two mirrors. One 
cycle is therefore equivalent to a two-lens system with two identical lenses and 
the same distance between the lenses. Let us rewrite the matrix corresponding 
to the one-lens system given in Eq.J2IJ. 

C=( X ~ 1 (48) 
V x x — 1 J 

Then one complete cycle consists of C 2 . For N cycles, the expression should be 
C 2N . However this calculation, using the above expression, will not lead to a 
manageable form. However, we can resort to the expressions of Eci. (|35J) and 
Eq.JsHJ). Then one cycle consists of 

C 2 = [B(ri)RW2)B(-ri)][B(r,)R(<l,/2)B(-Ti)] 

= B(r))R(ci>)B(-r)), (49) 
if the upper-right element is negative. If it is positive, the expression should be 
C 2 = [B(ri)X(x/2)B{-v)][B(Tj)X(x/2)B{-Ti)] 

= B( V )X( X )B(- V ). (50) 
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If these expressions are repeated N times, 

C 2N = B(r))R(N<t,)B(-r)), (51) 

if the upper-right element is negative. It it is positive, the expression should be 

C 2N = B( V )X(N X )B(-r 1 ). (52) 

As N becomes large, cosh(iVx) and sinh(iVx) become very large, the beam 
deviates from the laser cavity. Thus, we have to restrict ourselves to the case 
given in Eq. (|5TJl . 

The core of the expression of Eq. (|51|l is the rotation matrix 

R(Ncf>) = [R((j))] N . (53) 

This means that one complete cycle in the cavity corresponds to the rotation 
matrix R((f>). The rotation continues as the beam continues to repeat the cycle. 

Let us go back to Eq. 149|) . The expression corresponds to a Lorentz boosted 
rotation, or bringing a moving particle to its rest frame, rotate, and boost 
back to the original momentum. The rotation associated with the momentum- 
preserving transformation is called the Wigner's little-group rotation, which is 
related to the Wigner rotation commonly mentioned in the literature |l(Jj . 



8 Multilayer Optics 

The most efficient way to study multilayer optics is to use the S-matrix formal- 
ism We consider a system of two different optical layers. For convenience, 
we start from the boundary from medium 2 to medium 1. We can write the 
boundary matrix as |31| 

B(n) - ( cosh ( 7 ?/ 2 ) sinh(77/2) \ 
[V1 ^sinh(77/2) cosh(77/2) J ' [ ' 

taking into account both the transmission and reflection of the beam. As the 
beam goes through the medium 1, the beam undergoes the phase shift repre- 
sented by the matrix 

P(M={ o J l/2 )- (55) 
When the wave hits the surface of the second medium, the boundary matrix is 

B(-n)-( cosh (^/2) -sinh(ry/2)\ 
1 V) ^-sinh(77/2) cosh(77/2) J' [ ' 

which is the inverse of the matrix given in Eq.JMJ). Within the second medium, 
we write the phase-shift matrix as 

Pifo)=[ o e */ 2 j- (57) 
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Then, when the wave hits the first medium from the second, we have to go back 
to Ea. (|54ll . Thus, the one cycle consists of 

fcosh(7?/2) sinh(7 ? /2)\ / e -*/ 2 
\smh(r]/2) cosh^)/ I e^ 1 ' 2 



cosh(r?/2) -sinh(?7/2) \ f e^^ 2 
-sinh(r?/2) cosh(77/2) J \ e 1 ^/ 2 



(58) 



The above matrices contain complex numbers. However, it is possible to 
transform simultaneously 

/cosh( ?7 /2) sinhfo/2)\ 

^sinh(r?/2) coshfa/2),/ [ ' 

(60) 



to 



exp(?7/2) 
exp(-??/2) 



and transform 



to 



e -#i/2 o 

e 1 ^/ 2 



(61) 
(62) 



cos(&/2) - sin(&/2) 
sin(c/)i/2) cos(0i/2) 

using a conjugate transformation. It is also possible to transform these ex- 
pressions back to their original forms. This transformation property has been 
discussed in detail in Ref. [T4*| . 

As a consequence, the matrix of Ea. (|58|l becomes 

cos(0i/2) — sin((^i/2) 
e-i/ 2 J ^sin(^i/2) cos(0i/2) 

y (e-*>/ 2 \ / cos(0 2 /2) -sin(«fe/2)\ 

X ^ e^ 2 ^shu> 2 /2) cos(0 2 /2) J ' l0dj 

In the above expression, the first three matrices are of the same mathematical 
form as that of core matrix for the one-lens system given in Ea. (|35|l . The 
fourth matrix is an additional rotation matrix. This makes the mathematics of 
repetition more complicated, but this has been done |15| . 
As a consequence the net result becomes 

B(n)R(Na)B(-fi), (64) 

or 

BMXiNQBi-n), (65) 

where the parameters /i, a and £ are to be determined from the input parameters 
r), (f>\ and <\>%. Detailed calculations are given in Ref. |15| . 

It is interesting to note that the Lorentz group can serve as a computational 
device also in multilayer optics. 
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9 Concluding Remarks 



We have seen in this report that the Lorentz group provides convenient calcu- 
lational tools in many branches of ray optics. The reason is that ray optics is 
largely based on two-by-two matrices. These matrices also constitute the group 
SL(2, c) which serves as the universal covering group of the Lorentz group. 

The optical instruments discussed in this report arc the fundamental com- 
ponents in optical circuits. In the world of electronics, electric circuits form the 
fabric of the system. In the future high-technology world, optical components 
will hold the key to technological advances. Indeed, the Lorentz group is the 
fundamental language for the new world. 

It is by now well known that the Lorentz group is the basic language for 
quantum optics. Coherent and squeezed states are representations of the Lorentz 
group. It is challenging to see how the Lorentz nature of the above-mentioned 
optical components will manifest itself in quantum world. 

The Lorentz group was introduced to physics by Einstein and Wigncr to 
understand the space-time symmetries of relativistic particles and the covariant 
world of electromagnetic fields. It is gratifying to note that the Lorentz group 
can serve as the language common both to particle physics and optical sciences. 



Appendix 

A Lorentz Transformations 

Let us consider the space-time coordinates (t, z, x, y). Then the rotation around 
the z axis is performed by the four-by-four matrix 
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cos 9 — sin 9 

, sin 9 cos 9 



(66) 



This transformation is generated by 




(67) 



Likewise, we can write down the generators of rotations J\ and J2 around the 
x and y axes respectively. 



J 1 = 



0^ 

z 



,0 -i 0, 




(68) 
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These three generators satisfy the closed set of commutations relations 



(69) 



This set of commutation relations is for the three-dimensional rotation group. 
The Lorentz boost along the z axis takes the form 



/ cosh 77 
sinh?7 


V 



sinh r\ x 

cosh 77 

1 

1. 



(70) 



which is generated by 



fO i 0^ 

i 



Vo o, 



(71) 



Likewise, we can write generators of boosts K\ and K% along the x and y axes 
respectively, and they take the form 



K x = 



/0 i 0' 



1 

Vo 0, 



Ko = 



f0 i' 





1 0, 



These boost generators satisfy the commutation relations 



Km. 



[K, Kj 



(72) 



(73) 



Indeed, the three rotation generators and the three boost generators satisfy 
the closed set of commutation relations given in Ea. (|69l) and Ea. (|73|l . These 
three commutation relations form the starting point of the Lorentz group. The 
generators given in this Appendix are four-by-four matrices, but they are not the 
only set satisfying the commutation relations. We can construct also six two- 
by-two matrices satisfying the same set of commutation relations. The group of 
transformations constructed from these two-by-matrices is often called SL(2, c) 
or the two-dimensional representation of the Lorentz group. Throughout the 
present paper, we used the two-by-two transformation matrices constructed from 
the generators of the SL(2, c) group. 



B Spinors and Four- vectors in the Lorentz Group 

In Appendix A, we have noted that there are four-by-four and two-by-two rep- 
resentations of the Lorentz group. The four-by-four representation is applicable 
to covariant four-vectors, while the two-by-two transformation matrices are ap- 
plicable to two-component spinors which in the present case are Jones vectors. 
The question then is whether we can construct the four-vector from the spinors. 
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In the language of polarization optics, the question is whether it is possible to 
construct the coherency matrix |381 139j from the Jones vector. 

With this point in mind, let us start from the following form of the Pauli 
spin matrices. 

cti = (J -°i)< ct2 = (? i)> CT3= (* o*)- (74) 

These matrices are written in a different convention. Here 173 is imaginary, 
while (T2 is imaginary in the traditional notation. Also in this convention, we 
can construct three rotation generators 

Ji = (75) 
which satisfy the closed set of commutation relations 

[Ji, Jj] = ieijkJk- (76) 
We can also construct three boost generators 

Ki = (77) 

which satisfy the commutation relations 

[K i ,K j ] = -ie ijk J k . (78) 

The Ki matrices alone do not form a closed set of commutation relations, and 
the rotation generators Ji are needed to form a closed set: 

[J i ,K j ]=ie ijk K k . (79) 

The six matrices Ji and Ki form a closed set of commutation relations, 
and they are like the generators of the Lorentz group applicable to the (3 + 
l)-dimensional Minkowski space. The group generated by the above six matri- 
ces is called SL(2, c) consisting of all two- by- two complex matrices with unit 
determinant. 

In order to construct four-vectors, we need two different spinor represen- 
tations of the Lorentz group. Let us go to the commutation relations for the 
generators given in Eas. i|7t)[l . (|78|l and l|79|l . These commutators are not invari- 
ant under the sign change of the rotation generators J, , but are invariant under 
the sign change of the squeeze operators Ki. Thus, to each spinor representa- 
tion, there is another representation with the squeeze generators with opposite 
sign. This allows us to construct another representation with the generators: 

Ji = ^ = (80) 
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We call this representation the "dotted" representation. If we write the trans- 
formation matrix L of Eq. (|15f) in terms of the generators as 

^cri) j , (81) 

then the transformation matrix in the dotted representation becomes 

L = eJ-iT i (6j(Xi - WfcQj) j ■ (82) 

In both of the above matrices, Hermitian conjugation changes the direction of 
rotation. However, it does not change the direction of boosts. We can achieve 
this only by interchanging L to L, and we shall call this the "dot" conjugation. 

Likewise, there are two different set of spinors. Let us use u and v for the 
up and down spinors for "undotted" representation. Then u and v are for the 
dotted representation. The four-vectors are then constructed as |4L)j 

uii = —(x — iy), v v — (x + iy), 

uv = (t + z), vu = —(t — z) (83) 

leading to the matrix 

„ / uv —uu\ ( u\ , . . \ , a ., 
C=[ . • = (« ~u), (84) 
y vv —vu J \v J 

where u and ii are one if the spin is up, and are zero if the spin is down, while v 
and v are zero and one for the spin-up and spin-down cases. The transformation 
matrix applicable to the column vector in the above expression is the two-by-two 
matrix given in Eq. l|15|l . What is then the transformation matrix applicable to 
the row vector (v, —it) from the right-hand side? It is the transpose of the 
matrix applicable to the column vector (i), — ii). We can obtain this column 
vector from 

v 



(85) 



by applying to it the matrix 



9 = -iv 3 =(° 1 Q 1 ). (86) 

This matrix also has the property 

ga.g- 1 = - (a t ) T , (87) 

where the superscript T means the transpose of the matrix. The transformation 
matrix applicable to the column vector of Eg. (|85|) is L of Ea. (|82|l . Thus the 
matrix applicable to the row vector (v, — u) in Ea. (|84|l is 



{g-'Lg} 



T = g- l L T g. (88) 
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This is precisely the Hermitian conjugate of L. 

In optics, this two-by-two matrix form appears as the coherency matrix, and 
it takes the form 



r - ( < E i Ex > < E y Ex > | (qq\ 



where < E*Ej > is the time average of E*Ej. This matrix is convenient when 
we deal with light waves whose two transverse components are only partially 
coherent. In terms of the complex parameter w, the coherency matrix is pro- 
portional to 

C ={re~- ( 9 °) 

if the x and y components are perfectly coherent with the phase difference of 
S. If they are totally incoherent, the off-diagonal elements vanish in the above 
matrix. 

Let us now consider its transformation properties. As was noted by Opatrny 
and Perina [2Zj , the matrix of Eq. (|89(1 is like 

c= (t + z x-iy\ (91) 
\x + iy t-z J 

where the set of variables (x,y,z,t) is transformed like a four-vector under 
Lorentz transformations. Furthermore, it is known that the Lorentz transfor- 
mation of this four-vector is achieved through the formula 

C = LCL\ (92) 

where the transformation matrix L is that of Ea. (|15|) . The construction of 
four- vectors from the two-component spinors is not a trivial task |4(JI I41j . The 
two-by-two representation of Ea. l|91|l requires one more step of complication. 



C Conjugate Transformations 

The core matrix of Eq. contains the chain of the matrices 

(e-* \ /cosh 7/ sinh 77 \ f \ 

\ J \si-ahr] coshrj J \ e*« J ' 1 ' 

The Lorentz group allows us to simplify this expression under certain conditions. 

For this purpose, we transform the above expression into a more convenient 
form, by taking the conjugate of each of the matrices with 

i 



Then dWC^ 1 leads to 



cos <j> — sin <p \ ( cosh r\ sinh r\ \ / cos £ — sin £ 
sin <f> cos <fi J \ sinh 77 cosh 77 J \ sin £ cos £ 



(95) 
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In this way, we have converted W of Eg. l|93|l into a real matrix, but it is not 
simple enough. 

Let us take another conjugate with 

C2 =^(-i 0- (96) 

Then the conjugate C^CiWCTf becomes 



cos 4> — sin (j) \ / e n \ ( cos £ — sin £ 
sin0 cos0 Mo e~ v J \sin£ cos£ 



(97) 



The combined effect of C2C1 is 

C = C 2 Ci = —= ( -m/4 -in/i ) i (98) 



1 ( e™/ 4 e W 4 



-i7r/4 pin I A. 



with 

"Til 

After multiplication, the matrix of Eq. (|97() will take the form 



c = — ( /4 X/4 ) ■ (99) 



y =(c £)> (100) 

where A, B, C, and D are real numbers. If B and C vanish, this matrix will be- 
come diagonal, and the problem will become too simple. If, on the other hand, 
only one of these two elements become zero, we will achieve a substantial mathe- 
matical simplification and will be encouraged to look for physical circumstances 
which will lead to this simplification. 

Let us summarize, we started in this section with the matrix representation 
W given in Eq. 1|93[1 . This form can be transformed into the V matrix of Eq. (|97|l 
through the conjugate transformation 

V = CWC~ 1 , (101) 

where C is given in Ea. (|98fl . Conversely, we can recover the W representation 

by 

W = C- X VC. (102) 

For calculational purposes, the V representation is much easier because we are 
dealing with real numbers. On the other hand, the W representation is of the 
form for the S-matrix we intend to compute. It is gratifying to see that they 
are equivalent. 

Let us go back to Ea. Q97f) and consider the case where the angles 4> and £ 
satisfy the following constraints. 

+ £ = 20, </>-£ = 7r/2, (103) 
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thus 

(t> = 9 + TT/A, £ = 0-7r/4. (104) 
Then in terms of 6, we can reduce the matrix of Eg. l|97l) to the form 



(cosh i]) cos(2#) sinh r\ — (cosh r\) sin(20) 

sinh rj + (cosh 77) sin(20) (cosh 77) cos(20) 



(105) 



Thus the matrix takes a surprisingly simple form if the parameters 9 and 77 
satisfy the constraint 

sinh 77 = (cosh 77) sin(20). (106) 



Then the matrix becomes 

1 
2 sinh rj 1 



(107) 



This aspect of the Lorentz group is known as the Iwasawa decomposition , 
and has been discussed in the optics literature E] . 

The matrices of the form is not so strange in optics. In para-axial lens optics, 
the translation and lens matrices are written as 

1 1). ^ (i :)• < io8) 

respectively. These matrices have the following interesting mathematical prop- 
erty 0, 

1 u 1+ u 2 \ (1Q9) 



i 


Ml \ 


i 


T ■( 




1 ) 







n o 1 

and 

' 1 o\ / 1 o\ _ / 1 

ui 1 / V u\ 1 / \ ui+u 2 1 



(110) 



We note that the multiplication is commutative, and the parameter becomes 
additive. These matrices convert multiplication into addition, as logarithmic 
functions do. 



D Euler versus Lie Representations 

In this paper, we restricted ourselves to the algebra of two-by-two matrices and 
avoided as much as possible group theoretical languages. In order to explain 
where those algebraic tricks came from, we give in this Appendix A group 
theoretical interpretation of what we did in this paper. 

The group SL(2, c) consists of two-by-two unimodular matrices whose ele- 
ments are complex. There are therefore six independent parameters, and thus 
six generators of the Lie algebra. This group is locally isomorphic to the six- 
parameter Lorentz group or 0(3, 1) applicable to the Minkowskian space of three 
space-like directions and one time-like direction. 

Like the Lorentz group, the SL(2, c) has a number of interesting subgroups. 
The subgroup most familiar to us is SU(2) which is locally isomorphic to 
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the three-dimensional rotation group. In addition, this group contains three 
subgroups which are locally isomorphic to the group 0(2, 1) applicable to the 
Minkowskian space of two space-like and one time-like dimensions. 

One of the subgroups of SL(2, c) is SL(2, r) consisting of matrices with real 
elements. This subgroup is also called the Sp(2) group which we used in this 
paper in order to carry out the Iwasawa decomposition. Another interesting 
subgroup is the one we used for computing the S matrix, which starts with the 
boundary matrix of Ea. l|5^|l and the phase-shift matrix of Eq. l(531) . This group 
is called SU(1, 1). The present paper exploits the isomorphism between Sp(2) 
and SU(1, 1). While the physical world is describable in terms of SU(1, 1), we 
carry out the Iwasawa decomposition in the Sp(2) regime. 

Indeed, the conjugate transformation of Ea. H93l) to Ea. (|95|l is from SU(1, 1) 
to Sp(2), while the transition from Ea. l|95|l to Ea. (|97() is within the Sp(2) group. 
Thus, the transition from Eq. (JHSJ to Ea. (|97|l is a conjugate transformation from 
the SU(1, 1) subgroup to the subgroup Sp(2) of SX(2,r). 

In this paper, we are concerned with the decomposition of the Sp(2) and 
SU(1, 1) matrices. Unlike the traditional approach to group theory which starts 
from the generators of the Lie algebra, we used in this paper an approach 
similar to what Goldstein did for the three-dimensional rotation group in terms 
of the Euler angles [2SJ - There are three-generators for the rotation group, but 
Goldstein starts with rotations around the z and x directions. Rotations around 
the y axis and the most general form for the rotation matrix can be constructed 
from repeated applications of those two starting matrices. Let us call this type 
of approach the "Euler construction." 

There are three basic advantages of this approach. First, the number of 
"starter" matrices is less than the number of generators. For example, we 
need only two starters for the three-parameter rotation group. In our case, 
we started with two matrices for the three-parameter group Sp{2) and also for 
SU (1,1). Second, each starter matrix takes a simple form and has its own 
physical interpretation. 

The third advantage can be stated in the following way. Repeated applica- 
tions of the starter matrices will lead to a very complicated expression. However, 
the complicated expression can decomposed into the minimum number of starter 
matrices. For example, this number is three for the three-dimensional rotation 
group. This number is also three for SU(2) and Sp(2). We call this the Euler 
decomposition. The present paper is based on both the Euler construction and 
the Euler decomposition. 

Among the several useful Euler decompositions, the Iwasawa decomposition 
plays an important role in the Lorentz group. We have seen in this paper what 
the decomposition does to the two-by-two matrices of Sp(2), but it has been 
an interesting subject since Iwasawa's first publication on this subject |28|. It 
is beyond the scope of this paper to present a historical review of the subject. 
However, we would like to point out that there are areas of physics where this 
important mathematical theorem was totally overlooked. 

For instance, in particle theory, Wigner's little groups dictate the internal 
space-time symmetries of massive and massless particles which are locally iso- 
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morphic to 0(3) and E{2) respectively [7j. The little group is the maximal 
subgroup of the Lorentz group whose transformations do not change the four- 
momentum of a given particle 35 . The _E(2)-like subgroup for massless particles 
is locally isomorphic to the subgroup of SL{2, c) which can be started from one 
of the matrices in Eq. (|108|) and the diagonal matrix of Eq. I|55|) . Thus there was 
an underlying Iwasawa decomposition while the the i?(2)-like subgroup was de- 
composed into rotation and boost matrices |36| . but the authors did not know 
this. One of those authors is one of the authors of the present paper. 

In optics, there are two-by-two matrices with one vanishing off-diagonal el- 
ement. It was generally known that this has something to do with the Iwasawa 
effect, but Simon and Mukunda [33] and Han et al. jjj started treating the 
Iwasawa decomposition as the main issue in their papers on polarized light. 

In para-axial lens optics, the matrices of the form given in Eq. I|108[) are the 
starters 021 j and repeated applications of those two starters will lead to the most 
general form of Sp(2) matrices. It had been a challenging problem since 1985 [4*2"] 
to write the most general two-by-two matrix in lens optics in terms the minimum 
number of those starter matrices. This problem has been solved recently |1(J| . 
and the central issue in the problem was the Iwasawa decomposition. 
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